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One important way to characterize the states having a gauge symmetry spontaneously broken over 
multibands should be to look at their collective excitation modes. We find that a three-band system 
has multiple Leggett modes with significantly different masses, which can be classified into different 
dynamical classes according to whether multiple inter-band Josephson currents add or cancel. This 
provides a way to dynamically characterize multiband superconductivity while the pairing symmetry 
is a static property. 



PACS numbers: 74.20.-z,03.75.Kk 

Introduction — Superconductivity and superfluidity, 
with their spontaneous broken gauge symmetry, harbor 
unexpected fascinations. Specifically, the seminal discov- 
ery of the iron-based superconductor 1 has kicked off re- 
newed interests in multiband superconductivity. There, 
the gauge-symmetry breaking involves multiple bands, 
so the ground state is expected to possess the features 
that can occur only in multiband systems, such as the 
s± wave pairing proposed for the iron-based supercon- 
ductor 2 . Even more interesting are excited states, es- 
pecially the collective excitations of condensate phases 
associated with the broken gauge symmetry- . For two- 
band superconductors, there is a classical 1966 work by 
Leggett, who has shown that a two-band superconductor 
accommodates a special collective excitation, now called 
Leggett's mode, emerging as an out-of-phase mode be- 
tween different superfluids 4 . The Leggett's mode was 
experimentally detected in a two-band superconductor, 
MgB 2 (Ref.—) via the Raman scattering—. 

However, little is known as to what happens to the 
collective modes when there are three or more bands. In 
fact, despite intensive studies on the gap function sym- 
metries in the ground states for various classes of super- 
conductors, dynamics of collective excitations in multi- 
band systems has yet to be systematically investigated. 
The question is becoming increasingly realistic, since the 
iron-based superconductors, for instance, have been re- 
vealed to have a five-band structure, where three bands 
contribute to both the Fermi surface and the gap func- 
tion—, as subsequently experimentally examined—. To 
study the dynamics of collective modes in multiband sys- 
tems is thus expected to be a new avenue together with 
unconventional pairing symmetry. 

The purpose of the present paper is to examine, in gen- 
eral, whether a novel class of Leggett's modes indeed ex- 
ists specifically in superconductors and superfluids that 
have three or more bands, which should shed light on 
dynamics of multiband condensates. Naively, one might 
expect that the two- and three-band cases may be simi- 



lar. Here, however, we shall show, on the basis of an ef- 
fective action for the phase fluctuations, that multiband 
superfluids having three or more bands are in fact unique 
in that there exist multiple Leggett's modes classified by 
a dynamical class introduced here. The distinction be- 
tween the dynamical classes come from the presence of 
multiple inter-band Josephson couplings, whose additive 
and subtractive combination to the effective action are, 
respectively, classified as class "even" and "odd." The 
two Leggett's modes in the class odd are predicted to 
have significantly different masses, which is testable and 
serves to characterize multiband superconductors. 

Formulation — Let us start with a model for multi- 
band superfluidity /superconductivity. We take the sim- 
plest possible Bardeen-Cooper-Schrieffer type Hamilto- 
nian density for an TV-band superfluid, assumed to be 
neutral here, is given as ti(r) = Hkmetic(r) + Hp a frm S (r) 

with Akinetic = J2iLl Sa=t,4.^!<T e i( - ^)^«7 wnere tne 

dispersion £j(fe) is assumed to be parabolic with H = 1, 

While "Hpairing = ~J2i,j=l 9ij^i^i^ji^jt- Here i>i<r 

(i = 1, 2, . . . , N) is the field operator for the ith band in 
a superconductor or ith atomic species in an atomic gas, 
and Q = (gij) the pairing matrix, taken to be positive- 
definite, and we drop its fe-dependence. While Q may 
come from (the short-range part of) the electron corre- 
lation, we ignore the Coulombic part of the interaction, 
since the mass of Leggett's modes is not affected by the 
Anderson-Higgs mechanism, although the group velocity 
of the mode reflects the Coulomb interaction in super- 
conductors as shown in Ref. 8 for N — 2. While impairing 
with fc-dependence dropped, as is done in Leggett's orig- 
inal treatment, is too simple to examine complicated 
band structures as in the iron-based superconductors, we 
should initially know the basic physics for collective ex- 
citations in multiband superfluidity /superconductivity in 
the simplest case. 

The grand partition function is given, with 
the imaginary-time functional integral method, 
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in terms of a set of auxiliary fields AW as 
Z = Z /n i =i x,A(l)x,A(i) * e ^ Soff - Here the effective 
action is S cB = tfdr J d 3 r £^ =1 A^*(Q-% A« - 

X)i!=i hiGo.i — X)i!=i lnG^ 1 , where we assume 
the inverse Q 1 is well-defined (i.e., det£ > 0), Z a 
the partition function for the non-interacting system, 
and /3 the inverse temperature. Gi is Green's function 
in the Nambu representation for the ?th-band, which 
satisfies G~[ — Gq](I — Go,iKi), where Go,, is the 

free-fermion Green's function, / the unit matrix, and 
Ki = eriReA w — zo^ImA^, with a a being the Pauli 
matrices. The gaps are coupled through the inverse 
of the pairing matrix (Q^ 1 ), and we shall see that the 
non-zero off-diagonal elements of Q~ x determine the 
Leggett's modes. 

With this effective action the static gap equation reads 



N 

E 



tanh(/3£y/2) 
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where E t = ^J^ 2 + \W 2 , w c a cut-off frequency, and Nj 
the density of states (DOS) of the jth fermion on the 
Fermi surface. Here, we assume that each of the A W 's is 
constant (an s-wave). We can then look at the phase (f^ 

in Aj = |Ai|e lp o . The gap functions on different bands 
can take either the same or opposite signs, i.e., the phase 
difference, (f>Q = <Pq ^ — <fo\ takes either or ±7r and 
should obviously satisfy Y^iLi %' * = 0mod27r with 

Vo^" 1 " 1 ^ = Vo 1 '' as depicted in FigQTa). 

Let us derive an effective action for the fluctua- 
tions in the superfluid phase to single out the collec- 
tive dynamics at zero temperature. We first decom- 

(i) I '\ (i) 

pose the phase ip + <p (l> , into the equilibrium ip Q 
(as obtained in the mean-field gap equation) and the 
phase fluctuation A' ■ Around the solution of Eq. (fTJ , 
we obtain an action, S e g — dr J drVJnterband + 
Eti E~=i Tr (G' 0ji Ki) m /m, where 

Vinterband = ^ E Vij^ij I 1 _ COs(</7 (j) - <p W )], (2) 

% = cos(^' j) + Kij), A y - = 4|(g- 1 ) lJ ||A,||A J |, (3) 

comprises a sum of Josephson-couplings between the 
phases of different superfluid gaps that represents the 
inter-band Josephson currents caused by the rela- 
tive phase fluctuations. We mention that the inter- 
band Josephson couplings are analogous quantiles with 
the standard Josephson couplings but derived in non- 
perturbative manner. In the above expression, the con- 

or ±7r is used, and K%j(= or tt) the 
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quantities are gauge-transformed, i.e., G' oi (x;x') = 



A 12 sin(A ) - 
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FIG. 1: Schematic diagrams for N = 3 for (a) the differences 
in the superfluid phases in equilibrium and (b) the inter- 
band Josephson coupling for the phase fluctuation A around 
the equilibrium (where arrows represent the phase difference 
associated with the inter-band Josephson current). 
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TABLE I: Classification of the signs of the coupling for N = 3. 



x — (r, r), with the Dyson equation transformed into 

Parity in the inter-band Josephson couplings — Let 

us begin with a trivial two-band (N = 2) case, which is 

relevant to the superconductivity in MgB 2 — . Equation 

© then reduces to Vinterband = (r7i 2 Ai 2 /2)[l - cos((p (2) - 

(p^ 1 ')] with a single inter-band Josephson coupling. The 

sign of — (Q~ 1 )i2(= 512/ dot Q) is equal to that of 312. 

(1 2) 

Hence K12 = when 512 > 0, for which 4>o = 
for VS n terband to give a stable gap solution according to 
Eq. ((3]). Similarly, K12 = 4>q — n when gi2 < 0. Thus 
we always have 7712 = 1, which implies that the sign of g\2 
is totally irrelevant to the spectrum of collective modes 
in the two-band case, as noted in Refs.— 

We now turn to the three-band case with N = 3. We 
have three kinds of the inter-band Josephson currents as 
schematically depicted in Fig.QJb). In this case the set of 
the signs of the couplings (7712, V23, f]3i) can be classified 
into: (i) all the signs positive, (ii) two positive, one nega- 
tive, (iii) one positive, two negative, and (iv) all negative. 
We note that, when two or three fte's are negative, the 
Hessian matrices associated with VSnterband always have 
negative eigenvalues at (tpW, <p™ , tp^) = (0,0,0)—, 
which implies that the solution of Eq. (JTJ) is not a sta- 
ble minimum of S c s- We can thus exclude cases (iii) and 
(iv) [Table IJ. 

Let us have a closer look at {^ci } an d { K ij}> from 
which {rjij} in Eq. ^ is determined. Note first that 
{k^} is governed by Q^ 1 rather than £/, where, e.g., 
-(£7 -1 )i2 = (512533 - 513523)/ detg for N = 3. If we 
introduce a shorthand, H = (K12, K23, K31), we can define 
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FIG. 2: (Color online) Dispersion relations for cjl+ (red 
line), Wl_ (blue) and ojo (green) for (a) the class even with 
the values of gij displayed for which A2/A1 = —1.00 and 
A3/A1 = -1.01, and (b) the class odd with A2/A1 = -0.73, 
A3/A1 = -0.38. We set u c /k B T c = 5000 and g u = ±. The 
DOS (Ni) and the Fermi velocity are assumed to be identical 
between the bands. 



an integer, 



Pk = «i2 + K23 + «3i (mod 27T) 



(4) 



Since fty = or 7r, we have either pg/jt = or 1, which 
defines the parity "even" and "odd" classes, respectively. 

For the phases $q = (0q ',$0 i^o 3 '^)' we ^ ave f° ur 
possible cases, 0o = (0,0,0), (0,7r,7r), (0,7T, — n), or 



(1,2) 







(0, — 7T, 7r) as seen in Fig. [If a), where we assume 0q 
without a loss of generality. The first case corresponds 
to an s-wave with no sign changes, while the others sign- 

-rr j c A 1 , 2 ) , J.(2,3) . ,(3,1) 

reversing s- waves. It we detme p^ g = <p +0q +0q 
(mod 2tt), we have always p^ Q = 0. 

In the class even, in which all ryy's are 1, we find, for 
all the above possibilities for k and cj>Q, that the class 
even occurs only when pz = 0. On the other hand, the 
class odd, where one of ry^'s is —1, should satisfy pft = 1. 
Thus the class for {ry^} are completely characterized by 
the parity of p% alone, as summarized in Table [I] We 
remark that Tr (G- Kl) m appearing in S e s does not de- 
pend on {r]ij}, which implies that the effective action is 
completely distinguished by the dynamical class charac- 
terized by the parity of p^. Such dynamical classification 
should always be applicable in the pairing-interaction pa- 
rameter space when N > 3. In other words, {/%•}, which 

does not have to coincide with {</>□'} modulo 2-k, is no 
longer a matter of convention unlike the case of N = 2. 

Collective modes — We are now in position to calcu- 
late the collective modes. If we expand Seff around the 
stable solution of Eq. (QJ for y>W (g) with q = (ivi, q) and 
Vi = 2£tt/(3 (£ g Z) the Matsubara frequency, we have, 
to the leading order, S e s ~ ^Z q t{ Pq{M q / 4)0 q in the long- 
wavelength limit (|q| — > 0)—. Here we have introduced 
(p q = t ((p^(q),ip^(q), <p( 3 \q)), and a 3 x 3 real symmet- 



M 11 (q)+p 11 -T}i2\u ~V3ihi 

M„ = I -7yi2Al2 M 2 2(q) + fJ-22 -^23^23 

-%lA31 -??23A23 M 33 (q) + /J33 

where M u (q) = N^f + ^ W^PS.i/ 3 )?! witn u F£,i the 
£th component (£ = x, y, z) of the ith fermion's Fermi 
velocity. The diagonal elements involve the contributions 
from the inter-band Josephson currents, /in = ?yi2Ai2 + 
?73iA3i, etc. The dispersion relations for the collective 
excitation modes are given by the roots of det M q = 0. 
After an analytic continuation ivt — > co+iO, we obtain the 
three roots, wo, and wl-- Figure [2] displays the full 
dispersion. We can give explicit formulae for the three 
collective modes to the leading order in q for isotropic 
systems. The first root ujq corresponds to the Nambu- 
Goldstone mode with uj n = V\q\ + 0(|q| 2 ), where V 2 — 

E 3 =i(^Mot)«/3) and JVtot = ^. The other two 
roots, and are the Leggett's modes = 

4+4|9| 2 +0(|9| 4 ),where4 = ±U 2 Tv%V' 2 / {v 2 + -v 2 _). 
Here the "mass gap" (the Leggett's mode frequency at 
\q\ = 0, i.e., v+, vJ) are characterized by two quantities: 
The frequency scale is given by v 2 — r]i2V 2 2 + r]23 l/ ^ 
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^31^31 and vfj = (Ni+Nj)Xij/NiNj, where for Ay ^ 
corresponds to the frequency of the Leggett's mode in the 
two-band model. The difference v + — ^_(oc \/l — D) is 
characterized by 



D = 



4iVtot (^^^l A12A31 + ryi2?723Al2A23 + >723?73lA23A3l) 

D A NxN2N 3 



For D — > 1 the mass difference v+ — V- vanishes, whereas 
for D — > ojl- becomes a gapless mode (i.e., wl- — > 
c-\q\). We find that v + — v_ becomes large due to a 
cancellation among the terms in the numerator of D oc- 
curs for class odd. The term quadratic in \q\ for the two 
Leggett's modes are characterized by c±, which involve 
two velocities, U 2 - (1 - D)- 1 / 2 ELi(aW^ 2 )(*4/3) 
and V' 2 = v 2 /3 - V 2 , with v 2 = Ei=i u Fi- In thc 
limit where two bands are decoupled (A23 = A31 = 0), 
we recover the two-band result— with two gapless modes 
(u;o,u;l-), and one gapped mode along with v + 

and v- that coincide with those in Ref.— . 

Having formulated the Leggett's modes, let us see 
how they reflect the difference in the dynamical classes 
even and odd. If we look at the simplest case of iVj = 
Ni/N tot = i, and v 2 t = v 2 Jv 2 = \ (for 1 = 1,2,3), 
for gn = 322 = 333 = \ with gij = g^/lrQ, Fig[2] dis- 
plays typical dispersion relations for loq and wl± for the 
two dynamical classes. We immediately notice that the 
mass difference i>+ — V—(= wl+ — oj- at |q| = 0) is much 
greater in the class odd than in the class even. System- 
atic variation of v + — V- on an interaction parameter 
space (312,313) is shown, first in Fig.[3ja) for the case 
of two negative (repulsive) and one positive (attractive) 
inter-band couplings. The allowed parameter region is 
restricted by the positive-definiteness of Q—. The mass 
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FIG. 3: (Color online) Color-coded mass difference (oc 
VI — D) in the multiple Leggett's modes plotted on 312-531 
plane for Tja = i. We assume identical normal state proper- 
ties between the three bands, (a) Case of two negative (re- 
pulsive) and one positive (attractive) inter-band couplings. 
Yellow lines represent the boundaries between the even and 
odd classes, (b) Case of all the three inter-band couplings 
negative (repulsive), where only the class odd is allowed. In 
the blank area, the Leggett's modes become unstable. 



difference is close to zero in the central area of the class 
even region, while it is relatively large in the class odd 
region. Thus, we recognize that the mass difference of 
the Leggett's modes is an indicator of the difference in 
the dynamical classes. Note, however, that the distinc- 
tion becomes blurred when the three-band system ap- 
proaches a two-band behavior, which occurs when one of 
the off-diagonal elements of Q~ x is much greater than the 
others (e.g., around the corners of the class odd region 
with a triangular shape), so that there is no jumps at the 
boundary. 

Next, Fig.[3][b) displays the result for the case of all 
the three inter-band couplings negative (repulsive) with 
.912,313,523 < 0. Notably, this case always has the class 
odd throughout. This implies that the mass difference of 
the Leggett's modes is lower bounded. The result shows 
that the bound is ~ 0.7 for g 2 3 = —0.07 and gu = |. 
Even more interesting, we find that the Leggett's modes 



become unstable or ill-defined (within the present treat- 
ment which assumes real Aj's) in a narrow but finite 
region (blank area in Fig.O^b)) where all three (or two) 
of 512, 513, 523 take similar values. In this case, every two 
among the three gap functions Aj., A2, A3 in Eq. (QJ want 
to have opposite signs, but end up with complex values 
(i.e., the relative phase differences between the gaps de- 
viate from or tt) due to the "frustration", as pointed 
out in Ref.— . 

Summary — We have shown the presence of multi- 
ple dynamical classes in the N(> 3)-band superfluidity, 
which are characterized in terms of the parity of the mul- 
tiple inter-band Josephson couplings. We have revealed 
that the mass difference of the Leggett's modes is greater 
in the class odd, in contrast to the two-band case where 
the classification does not exist. So, the behavior of 
the multiple Leggett's modes is expected to characterize 
the dynamics of excitations in multiband superconduc- 
tors and superfluids. The expressions for the parities for 
N — 3 given here can be extended to N > 4, with their 
classification being an interesting future problem. Other 
future works include an extension of the present discus- 
sion to more general, fc-dependent case for realistic de- 
scription of multiband superconductors, such as the iron- 
based superconductors with their material dependence 2 , 
and an extension to anisotropic pairings. Another in- 
triguing problem is the Leggett's modes in the case where 
the gap functions break the time-reversal symmetry in 
"frustrated" three-band superconductivity 13 . Multiband 
superfluidity in cold fermionic atomic gases with multi- 
ple hyperfine states may also be an interesting playing 
ground where we have a greater tunability due to the 
Feshbach resonance^. How to detect Leggett's mode, for 
which Raman spectra have been discussed in Ref.— ~— , 
is also an important problem. 
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